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ABSTRACT 


A  finite  element  technique  has  been  used  to  determine  the  stress  intensity  at 
the  tip  of  a  bore  crack  in  a  cylindrical  fracture  specimen.  This  technical  note  describes 
the  detailed  input  for  the  computer  programs  which  perform  the  stress  analyses.  Much 
of  the  preparatory  work  involves  generating  a  suitable  mesh  of  elements  which 
represents  the  elastic  properties  of  the  specimen.  The  constraints  which  are  imposed  on 
the  mesh  geometry  are  described. 
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1.  INTRODUCTION 


Materials  Research  Laboratories  are  currently  investigating  the 
effect  of  residual  stresses  resulting  from  autofrettage  on  the  growth  of 
fatigue  cracks  originating  at  craze  cracks  on  the  bore  surface  of 
large-calibre  gun  barrels.  This  investigation  makes  use  of  ring 
specimens  25.4  mm  thick,  cut  from  76  mm  gun*  barrel  forgings;  the  rings 
are  compressed  across  a  diameter,  and  a  fatigue  crack  grows  from  a 
starter  notch  in  the  bore  surface  immediately  above  the  lower  loading 
point.  The  stress-intensity  at  the  tip  of  such  a  crack  is  known  for 
rings  having  a  diameter  ratio  of  0.5  [1,2);  in  this  case  however,  it 
was  necessary  to  use  specimens  with  a  diameter  ratio  of  0.555,  and  this 
factor  was  believed  to  have  a  significant  effect  on  the  stress  intensity 
calibration  of  the  cracked  specimen. 

To  obtain  the  stress-intensity  calibration  for  the  MRL  specimen 
geometry,  the  specimen  was  modelled  mathematically  and  subjected  to  a 
Finite-Element  (FE)  stress  analysis.  The  computer  programs  used  for 
this  work  were  made  available  by  ARL  and  are  known  collectively  as 
DISMAL  (D^placement  Method  anALysis) .  This  report  describes  the  use 
of  DISMAL  for  the  ring  specimen  geometry. 

2.  DISMAL 

The  DISMAL  programs  for  structural  analysis  were  developed  at 
the  University  of  Sydney  by  Dr  R.H.  Keays  (4).  Dr  Keays  later  worked 
at  ARL  and  modified  the  programs  for  operation  on  the  ARL  DEC  PDP-10. 
Facilities  for  inserting  a  special  crack  tip  element  were  adoed  later 
[3|  in  order  to  provide  better  representation  of  the  stress  gradient 
near  a  crack  tip.  Full  details  on  the  general  use  of  DISMAL  are  given 
in  reference  4. 
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2.1  Mesh  Generation 


To  analyse  a  specific  geometry,  it  is  first  necessary  to 
represent  the  structure  as  a  network  (or  mesh)  of  elements.  The  ring 
specimen  was  modelled  using  a  two-dimensional  mesh  consisting  of 
quadrilaterals  (QUAD4  elements),  triangles  (TRIM3)  and  one  special 
element  at  the  crack  tip.  The  symmetry  of  the  specimen  permitted  the 
use  of  a  mesh  representing  only  half  of  the  ring,  as  shown  in  Fig.  1. 

To  generate  the  mesh,  it  is  necessary  to  draw  out  the  specimen 
on  a  scale  sufficient  for  accurate  identification  of  the  element  corners 
(or  nodes),  including  features  such  as  the  fatigue  crack  and 
load-application  points.  The  special  element  is  then  defined  and 
consists  of  10  nodes  making  a  semicircle  centred  on  the  cr.  '.k  tip.  The 
node  points  are  separated  by  20“  and  are  placed  at  a  radius  'r'  where 
'r'  is  half  the  crack  length  or  half  the  remaining  ligament,  whichever 
is  smaller.  These  nodes  are  numbered  1  to  10  clockwise  and  are  linked 
consecutively.  Two  further  semicircles  of  10  nodes  are  positioned  at 
radii  l.lSr  and  1.30r,  and  serve  to  transmit  the  displacements  from  the 
body  of  the  structure  to  the  special  element;  these  are  numbered  11  to 
30  (clockwise,  increasing  radius).  These  nodes  are  linked  to  form  two 
bands  of  quadrilaterals  around  the  special  element. 

This  procedure  for  designating  the  corners  of  quadrilaterals 
(and  triangles  where  these  are  more  suitable)  is  continued  until  the 
whole  structure  is  represented.  The  semicircular  arrangement  of  points 
1  to  30  is  not,  in  general,  continued,  although  the  following  guide 
lines  are  followed. 

(a)  The  area  enclosed  by  each  element  should  gradually  increase 
•  th  distance  from  the  crack  tip  until  a  workable  element  size 
ti>,ically  0.5%  of  the  total  area)  is  achieved. 

(b)  The  numbering  system  used  should  ens.ire  that  the  difference 
between  the  largest  and  smallest  node  numbers  defining  an  element  is 
minimised  in  order  to  reduce  computing  time.  An  acceptable  band  width 
was  obtained  with  a  node  difference  of  15. 
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(c)  The  smallest  allovable  included  angle  is  40°  for  quadrilaterals 
and  15°  for  triangular  elements. 

(d)  Loads  should  be  introduced  via  a  series  of  elements  of 
gradually  Increasing  size. 

2 . 2  Input  to  Programs  -  File  "DATA" 

A  file  "DATA"  is  constructed,  containing  all  the  information 
required  to  characterize  the  mesh  and  loading  conditions.  The  format 
(including  order)  of  this  information  is  rigidly  defined  in  the  programs 
and  any  deviation  from  this  format  can  produce  results  which  contain 
errors.  Appendix  1  is  the  DATA  file  for  a  ring  specimen  which  contains 
a  crack  0.2  of  the  wall  thickness  in  length.  The  major  features 
requiring  care  during  construction  of  this  file  are  described  below. 

(a)  The  program  title  must  be  the  first  line  of  data;  a  FORTRAN 
format  of  18A4  is  used,  read  with  right  justification,  e.g.  RING 
_555_A/V0.2. 

(b)  Units  used  in  the  program  form  the  second  line,  in  2A3,  2X,  2A3 

format,  e.g.  MM _ NEWTON. 

(c)  MASK  is  the  third  line,  and  suppresses  the  unused  degrees  of 

freedom  for  each  node.  The  format  is  MASK  XYZ  R  R  R  where  XYZ 

-  —  -  X  y  z 

represents  displacement  in  the  x,  y  and  z  directions,  and  R^ ,  R^ ,  R^ 
allow  rotations  about  the  three  axes;  a  code  1  inhibits  movement  and  0 
allows  displacement,  e.g.  MASK  _  _  _  000111  allows  movement  of  the  nodes 
but  no  rotations. 

An  error  was  discovered  in  the  existing  program  prohibiting  the 
use  of  a  two-dimensional  mask  i.e.  001111,  and  it  was  therefore 
necessary  to  insert  the  'z'  direction  restriction  after  each  node 
number . 

(d)  The  fourth  line  defines  the  total  number  of  points  in  the 
structure  (and  hence  the  number  of  lines  following).  Format  is  POINTS, 
12X,  14. 
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(e)  The  following  lines  consist  of  information  about  each  node  in 
the  mesh;  number,  degree  of  freedom  and  position  (in  a  rectangular 
coordinate  system).  The  format  is  15,  2X,  611,  3F10.2,  and  nodes  are 
inserted  in  ascending  numerical  order  (in  fact,  the  program  re-numbers 
the  nodes  in  input  order,  so  it  is  necessary  to  take  care  with  input  to 
avoid  the  band-width  problem  outlined  earlier). 

In  the  example  in  Appendix  1,  it  should  be  noted  that  symmetry 
of  the  specimen  requires  the  use  of  a  y-direction  restriction  on  those 
points  along  the  symmetry  axis  and  common  to  both  halves  of  the 
specimen.  It  is  also  necessary  to  fix  one  point  on  the  mesh  in  space 
to  avoid  translational  movement  of  the  whole  specimen;  in  this  case, 
point  377  (on  the  symmetry  axis)  was  chosen. 

(f)  Following  the  last  node  description  are  two  lines  which  define 
the  material's  elastic  properties.  The  format  is 

MATERIAL,  17X,  TITLE  OR  COMMENT  (UP  TO  17  CHARACTERS) 

E15.0,  F5.3,  E15.0,  2E10.0 

where  the  second  line  contains  Young's  Modulus,  Poisson's  Ratio,  Shear 
modulus,  maximum  allowable  stress  in  tension  or  compression,  and  maximum 
allowable  shear  stress.  Default  for  the  maximum  allowable  stresses  is 
zero. 


(g)  The  next  section  defines  the  special  crack  tip  element;  two 
lines  are  used,  the  first  having  the  format  4A4,  A2,  IX,  14,  4A4,  A2 

SPECIALELEMENTS  _  1  _  CRACK_TIP 

The  second  line  is  a  mask  used  to  define  the  degrees  of  freedom  for  each 
node  point,  in  format  12,  6011,  the  first  term  being  the  number  of  nodes 
in  the  element,  the  others  being  masks  of  a  similar  form  to  those  used 
earlier.  Only  masks  of  000111  or  000000  can  be  used;  no  z-direction 
constraint  is  allowed.  The  third  line  (12,  1015)  defines  the 

similarity  and  corner  points  of  the  special  element;  a  second  special 
element  (identical  in  size  and  orientation)  may  be  inserted  by  making 
SIM=01  for  the  second  special  element.  Note  that  points  are  defined  in 
an  an t i -clockwise  direction. 
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(h)  Printing  of  the  element  stiffness  matrices  and  of  detailed 
stresses  may  be  switched  off  by  the  next  two  lines: 

PRINT_ELSTIF_OFF 

PRINTSTRESSOFF 

(i)  The  quadrilateral  element  control  line  and  element  definitions 
follow,  in  the  format 

4A4,  A2,  IX,  14,  2X,  4A4,  A2 
12,  415,  F10.4,  ElO.l,  F5.3,  ElO.l 

where  the  first  line  contains 

QUADRI _ QUAD4 _ 307 

The  second  4A4,  A2  is  available  for  comments  which  will  appear  in  the 
output  text,  and  the  second  line  contains  (in  order)  similarity,  the 
four  node  numbers  (anticlockwise),  thickness  (assumed  to  be  unchanged  if 
left  blank).  Young's  Modulus,  Poisson's  Ratio,  and  Shear  Modulus  (the 
last  three  defaulting  to  the  values  in  the  previous  MATERIAL  statement). 
The  second  format  is  repeated  for  each  element  of  the  specified  type. 

(j)  Triangular  element  control  line  and  definitions  follow  a 

similar  format  (with  315  stead  of  415).  Thickness  may  be  left  blank, 
defaulting  to  the  previously-defined  value. 

(k)  If  the  structure  has  been  fully  defined,  the  line  END 
STRUCTURE  is  inserted. 

(l)  Any  number  of  load  cases  may  now  be  inserted,  in  the  format 

LOADS,  14X,  14,  2X,  4A4,  A2 
15,  6F10.02 

where  the  first  line  contains  the  number  of  loads  in  the  load  case,  and 

the  load  case  name,  and  the  second  line(s)  contain  the  number  of  the 

node(s)  at  which  load  is  applied,  the  forces  in  x,  y,  z  directions,  and 
moments  M  ,  M  ,  M  . 

«  '  y  '  7. 
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(m)  The  line  END  _  _  LOADS  terminates  the  data  file. 

2 . 3  Input  to  Programs  -  File  "CDATA" 

A  second  input  file  is  required  to  define  the  properties  of  the 
crack  tip  special  element.  The  file  contains 

Number  of  nodes  in  the  special  element 
Young's  modulus,  Poisson's  ratio,  thickness 
and  radius  of  the  element 
Angular  separate  of  nodes  (in  radians) 

Integration  step  siae 

The  number  of  terms  in  the  stress  function 
Global  coordinates  of  local  axes  X.Y^Z,  (i=l,2,3). 

Further  details  of  these  parameters  are  given  in  reference  5,  and  an 
example  of  file  CDATA  is  given  in  Appendix  2. 

3.  DATA  CHECKING  AND  INPUT 

Two  programs  are  used  to  check  the  input  data.  The  first, 
DIPLOT  creates  a  plot  of  the  mesh  as  described  in  DATA,  and  this  plot 
may  be  inspected  for  continuity  and  geometric  anomalies.  An  example  of 
the  dialogue  required  to  run  DIPLOT  is  given  in  reference  4. 

The  second  program  INPUTS  performs  some  checking  in  preparation 
for  determination  of  the  stiffness  matrix.  If  there  are  no  errors  a 
communication  file  DTALK.TMP  is  created  on  disk.  If  no  such  file  is 
produced,  information  on  the  errors  in  DATA  is  given  in  file  PRINTl.LST. 

If  all  is  in  order,  the  programs  which  make  up  DISMAL  are  run 
as  a  BATCH  job  (too  much  core  is  required  for  running  the  programs 
individually  from  MONITOR).  The  batch  control  program  for  execution  of 
DISMAL  is  contained  in  Appendix  3. 
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4.  RESULTS  AND  DISCUSSION 


The  crack  tip  stress  intensity  is  given  in  Che  output  file 
PRINT3.LST.  For  this  particular  geometry,  the  symmetry  condition  used 
to  reduce  the  number  of  elements  implies  that  the  load  applied  to  the 
half-ring  is  half  of  that  which  would  be  applied  to  the  full  ring. 

Meshes  representing  crack  lengths  between  0.1  and  0.9  of  the 
ring  wall  thickness  (V)  were  processed.  The  values  of  stress  intensity 
(K)  obtained  for  a  load  P  are  given  in  Table  I,  in  the  form  of  the 
dimensionless  parameter  Y  where  B  is  the  ring  thickness  and 

It 

=  bTO 


For  ease  of  use,  a  polynomial  was  fitted  to  this  data,  giving 
the  relation 

Y  =  9.8054(a/U)'’ ■  ^  -  25.0044(a/W)  +  54.63f>6(a/U)^ '  ^  -  7B.330(a/V)^ 

*  62.0963(a/W)^ • ^  -  21 . 7418(a/W)’ 

This  polynomial  was  used  (6)  to  determine  the  relationship  between 
fatigue  crack  growth  rate  and  crack  tip  stress  intensity  in  ring 
specimens. 

5.  CONCLUSION 

A  finite  element  technique  has  been  used  to  determine  the 
stress  intensity  at  the  tip  of  a  bore  crack  in  a  ring  specimen  loaded  in 
compression  across  the  diameter;  various  crack  lengths  were  considered, 
and  the  results  have  been  presented  in  a  polynomial  form. 
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Table  I  Results  of  Analysis 


Crack  Length 

0.1  tf 
0.2  U 
0.3  W 
O.A  W 
0.5  U 
0.6  U 
0.7  V 
0.8  W 
0.9  U 


y  =  KB-Ju  /P 

1.717 

2.073 

2.273 

2.376 

2.433 

2.400 

2.312 

2.150 

1.865 
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Clack  tip  region  of  mesh  shown  in  Fig.  1.  Crack  tip  is  a 
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